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ABSTRACT 

The gamma-ray binary system PSR B1259-63/LS 2883 consists of a 48 ms pulsar 
orbiting a Be star. The system is particularly interesting because it is the only gamma- 
ray binary system where the nature of the compact object is known. The non-thermal 
radiation from the system is powered by the spin-down luminosity of the pulsar and 
the unpulsed radiation originates from the stand-off shock front which forms between 
the pulsar and stellar wind. The Be star/optical companion in the system produces 
an excess infrared flux from the associated circumstellar disc. This infrared excess 
provides an additional photon source for inverse Compton scattering. We discuss the 
effects of the IR excess near periastron, for anisotropic inverse Compton scattering 
and associated gamma-ray production. We determine the infrared excess from the 
circumstellar disc using a modified version of a curve of growth method, which takes 
into account the changing optical depth through the circumstellar disc during the 
orbit. The model is constrained using archive data and additional mid-IR observations 
obtained with the VLT during January 2011. The inverse Compton scattering rate was 
calculated for three orientations of the circumstellar disc. The predicted gamma-ray 
light curves show that the disc contribution is a maximum around periastron and not 
around the disc crossing epoch. This is a result of the disc being brightest near the 
stellar surface. Additional spectroscopic and near-infrared observations were obtained 
of the system and these are discussed in relation to the possibility of shock heating 
during the disc crossing epoch. 
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1 INTRODUCTION 

The gamma-ray binary system PSR B1259-63/LS 2883 con- 
sists of a 48 ms pulsar orbiting around a Be star. The 
orbit is eccentri c (e 0.87) with a n orb ital period of 
forb ~ 3.4 yr i|johnston et all Il992l . 1 19941) . The system 
is particularly interesting because it is the only gamma- 
ray binary system where the nature of the compact ob- 
ject is known. Unpulsed TeV gamma-ray emission has 
previously been detected around periast ron by the High 
Ener g y Ste reoscopic System (H.E.S.S.) l|Aharonian et all 
120051 . [20091 ) . Previous observations have also detected un- 
pulsed radio emission near periastron (e.g Ijohnston et al.l 
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Il999j ) and X-ray emission across the whole orbit (see 
e.g.lCominskv. Roberts fc Johnston1ll994l ; iKaspi et af]|l995l ; 



IChernvakova et al.l 12006 , 2009). The discovery of extended 
radio and X-ray str u ctures have been recently r e porte d 
<|Mold6n et all 1201 ll : iPavlov. Chang fc Kargaltsevl 1201 ll ). 
These observations have shown that the unpulsed emission 
is variable at all wavelengths with respect to the orbital pe- 
riod. The system was observed around the 2010 periastron 
passage (~ 14 December 2010) with Fermi/LAT in the high 
energy (HE, 0.1 — 100 GeV) regime and showed a flux mod- 
ulation during t he periastron passage, peaking ~ 30 days 
after periastron (|Abdo et al.ll201ll : rTam et al.ll201ll ). 

The non-thermal radiation from the system is powered 
by the spin-down luminosity of the pulsar and the unpulsed 
radiation originates from the stand-off shock front which 
forms between the pulsar and stellar wind. Extensive mod- 
elling of this process has been presented in lTavani fc Aronsl 
(| 19971 ). where the X-rays are produced via synchrotron radi- 
ation of ultra-relativistic electrons (Lorentz factor 7 ~ 10 6 ) 
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and gamma-rays through the inverse Compton (IC) scatter- 
ing of target photons from the Be star. The ratio between 
the X-ray and gamma-ray flux is dependent on the radia- 
tive and adiabatic cooling times of the post-shocked pulsar 
wind. IC cooling was predicted to increase near periastron. 
IC cooling of the pre-s hocked pulsa r wind was also suggested 
for PSR B1259-63 bv lBall fc Kirkl (|2000l V 



The optical companion in the system, LS 2883, is a Be 
star, and produces an infrared (IR) flux which is hig her than 
predi cted by Kurucz stellar atmosphere models (jKurucj 
1 19791 ) . This IR excess is produced by free-free and free- 
bound emission in the circumstellar disc which surrounds 
the Be star. These Be star circumstellar discs are decre- 
tion discs which form around Be stars and are regions of 
enhanced stellar outflow. Optical interferometry observa- 
tions of these extended circumstellar envelopes have shown 
that they are symmetric al around the rotation axis (see e.g. 
iQuirrenbach et al] [l994). The decretion discs are variable 
and grow and shr ink over periods of hundr eds to thousands 
of days. See e.g. iPorter fc Riviniusi ((2003) for a review of 
classical Be stars. 



The IR excess from the circumstellar disc provides an 
additiona l photon source to the blackbo dy flux from the op- 
tical star. IVan Soelen fc Meintiesl (|201ll ) showed that in the 
isotropic case this IR excess could increase the gamma-ray 
flux, particularly at GeV energies, by a factor ~ 2. It is pos- 
sible that near the disc crossing the IR excess may have a 
much greater effect. In this study we discuss the effects of 
the IR excess near periast ron, adopting an anisotropic in- 
verse Compton model (e.g. Fargion. Rostislav fc Salial 19971 ; 
iDubus. Cerutti fc Henrill2008l ). The soft photon excess pro- 
vides target photons for IC scattering in the Thomson limit 
where the scattering cross section is a maximum. The addi- 
tional effects of the differ ent electron cooling processes in the 
plerion shock front (e.g. iKirk. Ball fc Skiaeraasenlll999l). a s 
well as the effects of gamma-ray absorption ( |Dubusl l2006), 
will not be discussed in this paper and the reader is directed 
to the given reference for relevant discussions. In order to 
include the IR excess, we determine the target photon dis- 
tribution using a modified version of the curve of growth 
method which was first proposed by I Waters] (|l986l ) (referred 
to hereafter as the COG method) . This modified COG model 
takes into account the changing optical depth as seen by the 
pulsar during its orbit. We constrain the COG model using 
literature data and VLT observations which were taken dur- 
ing January 2011, approximately 22 days after the December 
2010 periastron passage. 



This paper is structured as follows: an outline of pre- 
vious observations and modelling of PSR B1259-63 is pre- 
sented in section [2] section [3] discusses observations we ob- 
tained with the VLT, SALT, SAAO 1.9m and IRSF tele- 
scopes after the 2010 periastron passage; section U dis- 
cusses the modified COG method that we have derived; the 
anisotropic IC scattering is discussed in section [5] section [6] 
outlines our modelling of the anisotropic inverse Compton 
scattering around periastron; sections [7] and [8] discuss the 
results and implications thereof. 



2 PSR B1259-63/LS 2883 

A few optical and IR photometric observations of LS 2883 
are reported in literature and available in bright star cata- 
logues. A comparison of the catalog ues' observation s show 
little change in the optical magnitude. I Johnston et ail (I1994T ) 
confined the spectral range of LS 2883 to a 09 - B2 type star 
with a colour excess of E(B — V) = 1.05 and assumed a spec- 
tral type B2e, with M* ~ 10 Mq and R, ~ 6 Rq. Analysis 
of spectroscopic observations also showed that H/3 emission 
occurred at a radius of 8.57?*, assuming a Keplerian circum- 
stellar disc, which im plies that the circums tellar disc has a 
radius i? d isc > 8.572, i|Johnston et al.lll994 ). 

More recently, iNegueruela et all (|20 1 if ) undertook spec- 
troscopic observations of LS 2883 and proposed new system 
parameters. These authors estimate an 09.5 V star with a 
mass of 31 Mq (for a non-rotating star), a radius of 9.2 Rq 
and an effective temperature of 33 500 K. The authors also 
suggest a new binary inclination of i = 33° and a lower 
colour excess of E(B — V) = 0.85. 

The size of the Be star's circumstellar disc can be con- 
strained by the observed pulsar eclipse around periastron. 
Observations of the system during consecutive periastron 
passages show that the pulsed radio signal is eclipsed be- 
tween roughly 20 days before periastron until 20 days af- 
ter. This is believed to be the result of the pulsar passing 
behind the Be star's circumstellar disc. The binary sepa- 
ration at ~ 20 days is ~ 45i?* (assuming M» = 31 Mq 
and i?* = 9.2 Rq) which implies that the disc radius is 
-Rdisc > 45.R* . For this reason we have adopted a disc radius 
of -Rdisc = 50-R* in this study. It should be noted, however, 
that Smooth Particl e Hy drodynamic (SPH) s imulations by 
lOkazaki et~ai1 i^OHl ) and iTakata et all (|2012f ) suggest that 
the circumstellar disc might be "blown-away" by the pul- 
sar wind if the base density of the disc is low. However, for 
larger densities, the disc remains more stable and these au- 
thors suggest that the density near the base of the disc is 
Pd ~ 10 -9 g cm -3 . This study only considers a constant disc 
size. 

PSR B1259-63 was detected by Fermi/L AT around 
the 2010 periastron with unexpected results. While the 
source remained faint before and just after periastron, ap- 
proximately 30 days after pe riastron, the sou r ce became 
very active and much brighter (|Abdo et al.ll201ll ;l Tarn et al.l 
2011). This was unexpected as it was assumed that the sys- 
tem would show fairly symmetric emission around perias- 
tron. Radio and X-ray observations around the same pe- 
riod showed no unusual or flaring event s and were consis - 
tent with previous periastron passages (|Abdo et alj|201lh . 
Significantly, the peak luminosity detected by Fermi/LAT 
was L-y fa 8 x 10 35 erg s _1 , which corresponds to approxi- 
mately 100 per cent of the spin-down power of the pulsar. No 
pulsed gamma-ray emission was detected with Fermi/LAT 
from PSR B1259-63 during this period. This is unusual as 
PSR B1259-63 fits most o f the criteria norm ally associated 
with gamma-ray pulsars <|Abdo et al.ll201ll ). The lack of a 
pulsed gamma-ray signal from a known radio pulsar poses 
interesting questions regar ding the location of gamma-ray 
pulse production. See e.g. lAbdo et al.l (|201p| ') for a discus- 
sion of the gamma-ray pulsars detected by Fermi/LAT. 

Possible explanations for the Fermi/LAT flare have 
been suggested in terms of IC scattering of the cold pul- 



Anisotropic inverse Compton scattering in PSR B1259-63 3 



sar wind (jKhangulvan et a or Doppler shifted syn- 

chrotron emission from the tail of the pulsar bow shock 
l|Kong et alJl2012h . 



3 OBSERVATIONS 

We obtained mid-IR observations of PSR B1259-63/LS 2883 
with the VLT telescope during January 2011. A number of 
observations were also undertaken at the SAACQ Suther- 
land observatory: we obtained near-IR observations with the 
IRSF telescope during December 2010, and spectroscopic 
observations with the SAAO 1.9m telescope and with SALI0 
in April 2011 and May 2012 respectively. The observations 
are summarised below. 



3.1 VLT observations 

We obtained (in Service Mode) mid-IR observations of 
PSR B1259-63/LS 2883 with the VL T Imager and Spec - 
trometer of the mid Infrared (VISIR; lLagage et all 12004 ') 
with the VLT on 5 January 2011. This is approximately 
three weeks after the 14 December 2010 periastron passage. 
The observations consisted of one pointing in multiple bands 
in the broad N-band region (8 - 13 /im). The observations 
were executed with a perpendicular chop/nod direction, a 
chopping amplitude of 10" and a random jitter of 3". The 
total integration time was of the order of 60 seconds per 
filter and achieved S/N ratios of < 10. The observations 
were calibrated with the standard star HD92682, which was 
observed as part of the standard ESO procedure, allowing 
for an average calibra tion accuracy within ~ 3 per cent. 
jSmette fc Vanzill2007( h 

The calibration star data were reduced with the 
visir_img_phot pipeline and the science target observa- 
tions were reduced using the standard visir_img_combine 
pipeline to correct for the image nodding and chop- 
ping. Both reduction pipelin es produce a single FITS 
l|Wells. Greisen fc Hartenlll98ll ) file with four measurements 
of the targets contained within each file. We calculated a flux 
conversion ratio for each filter by comparing the measured 
flux from the standard star to the defined standard flux. The 
flux conversion ratio was used to calibrate the measured flux 
from the science target. 

The final error in the VLT observation of LS 2883 was 
calculated by combining the estimated photometric error of 
each flux measurement and the standard deviation of the 
four measurements in each FITS file, for both the science 
target and standard star. Not included in the error calcula- 
tion is the ~ 3% error associated with the calibration process 
and any error in the standard star's catalogue flux value. 
The calibrated (but not de-reddened) flux at each filter is 
summarized in Table [1] 



1 South African Astronomical Observatory 

2 Southern African Large Telescope 



Table 1. VLT/ VISIR data for LS 2883 observed on 5 January 
2011. Flux measurements were calibrated using the standard star 
HD92682, but are not corrected of interstellar reddening. 



Filter 


central wavelength 


Flux 




(an 


mjy 


PAH1 


8.59 


261 ± 40 


Arlll 


8.99 


241 ± 30 


SIV 


10.49 


218 ± 32 


SIV.2 


10.77 


188 ± 21 


PAH2 


11.25 


208 ± 31 


SiC 


11.85 


198 ± 18 


PAH2.2 


11.88 


196 ± 54 


Nell 


12.81 


143 ± 51 



3.2 SAAO observations 

3.2.1 Spectroscopic observations 

We obtained spectroscopic observations of the Ha emission 
line from PSR B1259-63/LS 2883 on 10 April 2011 with 
the SAAO 1.9m telescope and on 6 May 2012 with SALlJfl 
(Fig. [TJ. The SAAO 1.9m observations were taken using a 
grating spectrograph (1200 lines mm -1 grating) which gives 
a ~ 1 A resolution. The data was reduced and analysed 
using the standard iraf/noao packages^ The SALT obser- 
vations were performed with the 0.6" slit and the pg2300 
grating, achieving a resolution of R ~ 10000. The SALT 
data was pre-reduced with the standard SALT pipeline 
|Crawford et alj|20ich and the data analysis was done with 
the standard iraf/noao packages. The Ha emission line is 
prominent in both observations and there appears to be no 
indication of stellar Ha absorption. This indicates that the 
circumstellar disc completely dominates the stellar absorp- 
tion at the same frequency. The observed equivalent width 
for the Ha line was EW = -50 ± 5 A in April 2011 and 
EW = -49 ± 2 A in May 2012. This indicates that the Be 
star was still surrounded by a substantial circumstellar disc 
after the disc crossing epoch. The observed equivalent width 
in April 2011 is in line with the measurement reported by 
iNegueruela et al.l (|201ll ) (EW = -54 ±2 A), while the emis- 
sion line is smaller in May 2012. An equivalent w i dth o f 
EW = 40 A was reported earlier by I Johnston et all (1 19941 ). 
and the smaller equivalent width detected in May 2012 most 
likely represents the normal variability associated with Be 
stars. 

3.2.2 IRSF 

Observations were undertaken with the InfraRed Survey 
Facility (IRSF) on 19 December 2010. The IRSF is a 
1.4 m near-IR telescope, located at SAAO near Sutherland 
and was built as a collaboration between SAAO and the 
Nagoya University in Japan. Using the SIRIUS3 camera, 
the IRSF takes simultaneous J, H and K s observations. The 

3 Program: 2012-1-RSA-003 

4 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under coop erative agree ment with 
the National Science Foundation. See e.g. iTodvl (1993). 

5 Simultaneous-3color InfraRed Imager for Unbiased Survey 
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6550 6555 6560 6565 6570 6575 6580 
Wavelength (A) 

Figure 1. Ha emission line of PSR B1259-63/LS 2883 taken 
10 April 2011 with the SAAO 1.9m telescope and on 5 May 2012 
with SALT. The wavelength range has been corrected to the he- 
liocentre. 



Table 2. IRSF observations obtained during December 2010 



MJD 


J 


H 


K s 


55549 


7.999 ± 0.038 


7.584 ± 0.045 


7.176 ±0.082 



data was reduced using the standard IRSF data reduction 
pipeline and data fro m the 2MASS point-source catalogue 
ijSkrutskie et alj|2006l ) was used to self-calibrated the frame. 
The observation data is summarized in Table [2] The final 
error is calculated by combining the photometric error es- 
timation, the measured error in the 2MASS data and the 
error in the self-calibration fit. The data is in accordance 
with the measurements obtained in the 2MASS catalogue. 



4 THE CURVE OF GROWTH METHOD 

IWatersI ([1986) proposed a COG method to model the IR ex- 
cess of Be stars, by determining the optical depth through 
the circumstellar disc. This model is appropriate for stars 
where the discs are viewed face-on at distances which are 
much larger than the size of the disc. Below we present a 
modified version of this COG method which we used to 
model the IR excess originating from the circumstellar disc 
as observed from the pulsar. 

In the COG model it is assumed that the circumstellar 
disc has a half-opening angle (?disc> extends to a radius J?disc, 
and has a power-law density profile that decreases with dis- 
tance as 



p(r) = po ) = /M/' 



where f = r/R, is the radial distance in units of the stellar 
radius (R*). It can be shown that the optical depth through 
the disc is given by 



r = iodise / r 2n d/, 



where i?„,disc is the o ptical depth p arameter for the circum- 
stellar disc given by |Waterslll986T ) 

-Ei/.disc = X\X,d- 

Here X\ is a wavelength dependent term, 

X x = A 2 {(1 - e- hv/kT ^)/{hu/kT disc )} 

x{g(v,T disc )+b(v,T di3C )} (2) 

where A is the wavelength, v is the frequency, k is the 
Boltzmann constant, Tdi sc is the disc temperature, and 
g(v, Tdiac) + b{v, Tdisc) is the sum of the free-free and free- 
bound gaunt factors. The variable which is fit by the 
COG method, is given by 



X* 



4.923 x 10 3 Vr d 7 s 3 c /2 AT 2 Wo 



R, 
Rq 



where z 2 is the mean of the squared atomic charge, /j, the 
mean atomic weight (in units of proton mass) and vj is the 
ratio of the number of electrons to the number of ions. 

We calculated the gaun t factors in equatio n ([2| follow- 
ing the method outlined in IWaters fc Lamersl (|l984h . This 
allows us to determine gaunt factors for a larger wavelengths 
range than was included in the published tables. 

The optical depth through the disc, as observed from 
any arbitrary point (along the direction ee,^ in Fig. [2]), can 
be defined in terms of an impact factor q, such that 



-2 ,2.2 

r = I + q , 



where all variables are in units of R, . Thus it follows that 



T = E V 



1//2 dr. 



For any value of q, the substitution w — arccos (q/r) can be 
made and it follows that 



r(q) = F^disc q 



-2n+l 



(3) 



where the integration is performed over suitable limits of 
w. The limits on the angle w are determined by the angle 
between the impact factor q and the positions where the line- 
of-sight intercepts the disc. An example is shown in Fig. [5] 
where the limits are given by the angles wi and W2- For a 
disc observed face-on at infinity q lies along the plane of the 
disc f or all lines-of- sight and the limits on w axe —9& lsc to 
"disc (|Waterdll986T) . For an observer close to the disc the 
value of r must be calculated differently for every line-of- 
sight as the optical depth will vary. The calculation of the 
limits of w and the value of q are presented in Appendix lAl 
It is assumed that the circumstellar disc is isothermal, 
and therefore the source function of the plasma in the disc 
can be described by Planck's function. The specific intensity 
of the radiation fiel d which originates from the disc is then 
determined by (e.g. IWright fc Barlow! Il975h 



(1) 



h = B„(Tdisc) (l-e" 



r(q)\ 



(4) 
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Figure 2. Line-of-sight through the circumstellar disc of a Be 
star. The centre of the pulsar /electron cloud is at O and the 
centre of the star is at position S. The optical depth along a line- 
of-sight through the disc is dependent on the length between the 
two disc intercepts xi,yi,zi and £2>?/2>Z2- See text for details. 



where -B„(Tdi ac ) is the Planck function, Tdi sc is the temper- 
ature of the disc and r(q) is determined by equation @. 

The unmodified COG model, as presented by IWaterj 

(1986) assumes that the circumstellar disc is viewed face- 
on, at a distance much greater than the size of the disc. The 
impact factor, q, then lies along the mid-plane of the disc 
and has a maximum value of q = Rdisc/R*- For q > 1 the 
specific intensity is determined by equation Q. For q ^ 1 it 
is assumed th e disc has a minimal influence on the intensity 
|WatersHl986h and then 

where I„ ;4 is determined from the appropriate Kurucz atmo- 
sphere or blackbody distribution. The ratio of the observed 
flux arisi ng from the disc to the flux arising from the star is 
given by (|Telting et al.lll998h 



" f.disc 



-Rdisc/-R. 



[1- 



,(9)1 



2q dq. 



(5) 



The equation above is appropriate to the flux which would 
be observed from Earth for a disc which is observed face-on. 
The COG method can then be used to model a Be star's 
optical/IR flux by modifying a Kurucz model atmosphere 
at IR wavelengths with the flux ratio given by equation ([5]) , 
which is calculated by fitting the five parameters n, X„, 
-Rdis c, Tdisc a nd fld isc to IR data. 

IWaterj (|l986h showed that the observed flux from a 
moderately inclined disc, viewed from a very large distance, 
will not deviate greatly from that of a disc observed face-on. 
For this reason we have calculated the parameters for the 
COG method using the unmodified model. However, when 
the pulsar is near the circumstellar disc, the disc target pho- 
ton spectrum, as observed from the pulsar's position, cannot 
be calculated assuming the disc is viewed at infinity. The 
optical depth, and hence the specific intensity and photon 
spectrum, must be calculated for all lines-of-sight through 
the disc from the position of the pulsar. This will lead to a 
variability in the observed photon spectrum as the optical 
depth changes during the orbit. 

In the inverse Compton calculation we have calculated 
the photon density arising from the circumstellar disc by 
using the modified COG method outlined above. 



5 INVERSE COMPTON SCATTERING 
5.1 Anisotropic scattering 

In order to model the effect of the circumstellar disc on IC 
gamma-ray emission in Be-XPBs, it is necessary to take into 
account the IR excess, the solid angle over the disc and the 
changing scattering angle. This is achieved by calculating 
the anisotropic IC at all points in the orbit. 

Anisotrop i c scat tering has been considered in detail by 
iFargion et al.l d 1997T) and has been adapted to the case of 
LS 5039 bv lDubus et ail l|200ct ). The total scattering rate is 
given by 



dJVtot 
dtdei 



«. P h(eo)n e (7) 



di\r 



SI J~f Je 



dt dei 



cosi?deo d7dS7, (6) 



where n p h and n are the photon and electron number dis- 
tribution respectively, diV 7 , eo /dtdei is the anisotropic scat- 
tering rate and the cosi? term corrects for the projected 
surface area of the circumstellar disc. The full equation for 
the Klein-Nishina scattering rate, dN lieo /dtdei, is derived 



and given in iDubus et al.l (|2008l . equation A6 in their pa- 
per). The total scattering rate is calculated by integrating 
over the appropriate range of electron and photon energies 
and the solid angle subtending the target photon source. 

5.2 Integration over the solid angle 

The solid angle subtending the target photons (star and/or 
disc), from the pulsar's perspective, changes during the or- 
bit. This, combined with the change in the scattering angle, 
results in a modulation of the IC scattering rate. For the case 
of the photon contribution from Be stars, the contribution 
from the disc must also be considered and the integration 
must occur over the whole disc. This is accomplished by 
summing the scattering contribution from the star and the 
disc. 

The solid angle subtending the star is given by 



dfi = sin 9 ( 
where <j> £ [0, 2ir], € [0, a*], and 

a* = arcsin I — J . 

The size of the solid angle is then dependent on the binary 
separation, r, between the star and the pulsar. 

In order to calculate the solid angle subtending the disc, 
the co-ordinate system shown in Fig. [5] is used. (See e.g . 
IPomme et ai1l2003l ; iMoskalenko fc Strondlioool ; |Trvkalll997j . 
for a general discussion of the formulae below.) 

The pulsar is at O which is at a vertical height H above 
the disc mid-plane and shifted a horizontal distance L (along 
the a;-axis) from the centre of the disc/star. The solid angle 
subtending the disc from the position O can be calculated 
by dividing the integration into two regions, Qi and Q2 '■ 



fir = / / sin(9d69dc/>, (L < R diB 
Jo Jo 

= 0, (L > R disc ), 



(7) 



where 



?i = arctan 



J?disc — L 



H 
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Figure 3. Coordinate system used to integrate over the solid 
angle subtended by a disc. 



and 

/ sinfld^dc?, 



(«) 



where 



62 = arctan 



/ itdisc + L \ 

V H ) 



and 

<P4>{9) — arccos 



Rj isc — L 2 — H 2 tan 2 
2Lif tan<9 



The total solid angle of the disc is = fii + Q2- 

The above equations apply to the solid angle over a flat 
disc (©disc = 0°), but can easily be extended to a circum- 
stellar disc with a finite half-opening angle #disc- The solid 
angle is given by equations (JT)) & (J8]) where the new limits 
on the polar and azimuthal angles are 



= tan 



tan 



1 / L — itdisc \ 
\ H — /tdisc / 



MO) 



H - h di , 

-Rdisc + L 
\H — /idisc 

R 2 isc -L 2 -(#-fe dis c) 2 tan 2 i 
2L{\H - /idi sc |)tanr? 



Here the additional term, /idisc = fldisc tan 6 disc , accounts 
for the opening angle of the disc. 



5.3 Scattering angle 

The anisotropic IC (AIC) scattering rate is calculated by 
assuming that the electrons are highly relativistic and that 
the photons scatter in the direction of the electrons' original 
momentum. This assumption was in corporated into the AIC 
scattering kernel (jDubus et al . 2008), and the scattering rate 
is dependent on the angle between the unperturbed path of 
the photon and the direction of the un-scattered electron, 
in the observer's reference frame. This angle, defined here 
as the scattering angle, 8 SC , is given by the angle between 
the direction of the unperturbed photons path (ee,^) and 
the direction of the observer (e bs), and is determined from 
(Fig.© 

e p h • e obs = cos# ac . 



Table 3. Combined VLT and archive data for LS 2883 de- 
reddened using E(B-V)=0.85. 



Band 


Wavelength 


Fo 




/ill! 


10™ 24 erg/s/cm 2 /Hz 




U.ODO 


CCfl HQ _|_ n 1 07 

ou.uo nz u. iu { 


LJ 


n aa 


^9 A'^'X -I- fl D££ 
DZ.^oo HZ U.UOo 


V 


0.55 


40.539 ± 0.097 


J 


1.235 


19.614 ±0.066 


H 


1.662 


12.695 ± 0.091 


K 


2.159 


10.677 ±0.029 


A 


8.23 


2.71 ± 0.14 


C 


12.13 


10.91 ± 5.88 


PAH1 


8.59 


2.64 ±0.41 


Arlll 


8.99 


2.43 ± 0.30 


SIV 


10.49 


2.19 ± 0.32 


SIV.2 


10.77 


1.89 ± 0.22 


PAH2 


11.25 


2.09 ± 0.31 


SiC 


11.85 


1.99 ±0.18 


PAH2.2 


11.88 


1.97 ±0.54 


Nell 


12.81 


1.43 ±0.52 



6 MODELLING THE AIC SCATTERING 

To calculate the IC scattering rate the infrared excess from 
the circumstellar disc must first be determined. 



6.1 Modelling the infrared excess 

We fitted the un-modified COG model to available archive 
data and new VLT observations. Archive optical data is 
available from a catalogue of bright stars near the Souther n 
Coalsack region compiled by IWesterlund fc~G arnicr (1989), 
while IR data is av ailable from the 2 MASS ( Skrutskie et al.l 
l200fj) and MffAfl jprice et all 120011 ) missions, as presented 
in I Van Soelen &: Meintiesl l|201lTl . The combined VLT and 
archive data for LS 2883 were de-reddened with the DISPO 
software pac kages using the co l our e xcess, E(B-V) = 0.85, 
proposed by iNegueruela et alj (|201ll ). The total combined 
and de-reddened data are shown in Table [3] 

Spectroscopic ob servations undertaken by 
INegueruela et all ([201 1 ) suggested a stellar tempera- 
ture of T* = 33 500 K and an effective gravity logg = 4.0 
for LS 2883 if the star was non-rotating. For this reason 
we fitted a Kurucz stellar atmosphere with a temperature 
T„ = 33 000 K and an effective gravity log g — 4.0 to the 
optical data points. The Kurucz model was fitted to the 
optical data at frequencies above 10 13 Hz to minimize the 
influence of the IR excess on the fit. 

Recent observations of Be stars suggest that the 
circumstellar discs are very thin (e.g. lOuirrenbach et al.l 



1 19971; IWood. Biorkman fc Biorkmanl Il997l; ICarciofil |2011mT 



In the SPH modelling presented by Okazaki et alj ( 201 lh 
& iTakata et al.l \20lj ) the disc is assumed to have a tem- 
perature of Tdisc = 0.6 T, and the height of the disc is as- 
sumed to change with radius as H(R,)/R* = 0.024 (imply- 
ing #disc ~ 0.7°). Therefore, we adopted a disc temperature 
of Tdisc = 0.6T* and a disc half-opening angle of 1°. 

6 Data from 2MASS and the Midcourse Space Experiment 
(MSX) are available on-line through the NASA/IPAC Infrared 
Science Archive. 
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Figure 4. Kurucz atmosphere (T* = 33 000 K, logg = 4.0) and 
COG method fitted to observational data of LS 2883 (Table 0. 
See text for details. 



We calculated the Gaunt factors in equation ([2|) assum- 
ing a plasma composed of 70 per cent hydrogen (H), 27 per 
cent helium (He) and 3 per cent carbon, nitrogen and oxygen 
(CNO). We assumed that the H and CNO is in an ionized 
state while the He is neutral. This was done in order to 
be con sistent with the data tabulated in lWaters fc Lamersl 
{ 19841 ). 

The COG model was fitted to the observ ational data 
using the Levenberg-Marquardt method (e.g. IPress et al.l 
Il98d ). Only the parameters n and X* were fitted while Tdi sc , 
#disc and .Rdisc were kept at assumed values. The parameters 
were assumed to have the values Tdisc = 0.6 T* = 19 800 K, 
#disc = 1° and i?di sc = 50-R„. 

Two COG fits were considered, one for free-free+free- 
bound scattering and the other for only free-free scattering. 
Fig. [3] shows the stellar atmosphere (dotted line), the pre- 
dicted IR flux from the disc for free-free (thin solid line) 
and free-free+free-bound (thin dashed line) scattering, and 
the total expected flux fitted to the data points in Table 
(thick solid and dashed lines). The flux from the MSX Band 
C (open square), which lies far above the predicted flux, was 
not included in the calculation of the COG fit as it appears 
to be significantly inconsistent with other measurements at 
this frequencies. The reported quality for Band C is 1, indi- 
cating a low significance. The Band C measurement there- 
fore represents an upper limit and not a true measurement. 
The detection in Band A (8.2 (im) is consistent with the ob- 
servations performed with the VLT, and is included in the 
fit. The best fit parameters found for the COG method are 
given in Table 3] 

The calculated free-free+free-bound Gaunt factors cre- 
ate large discontinuities in the predicted flux at frequencies 
> 2 x 10 14 Hz. At frequencies below this, the IR spectrum 
is dominated by free-free scattering. These large discontinu- 
ities are not s een in COG fits to other Be stars discussed in 
IWatersI (|l986l ). as these authors calculated the COG model 
to a maximum frequency of ~ 3 x 10 14 Hz, which is the 
limit of the tabulated Gaunt factors in IWaters fc Lamersl 
l|l984l ) . Below this frequency the influence of the free-bound 
scattering is still minimal. The free-free scattering pro- 
duces a smoother curve and this spectrum was adopted 



Table 4. Parameters of the COG method shown in Fig. [4] for 
free-free and free-free+free-bound (ff+fb) Gaunt factors, using 
the observational data in Table [3] 



Parameter 


free-free 


free-free+free-bound 


T* 


33 000 K 


33 000 K 


logs 


4.0 


4.0 


n 


3.0550 


2.8389 


log 10 X* 


10.245 


9.9073 


^disc 


50 R* 


50 R t 


-^disc 


19 800 K 


19 800 K 


#disc 


1° 


1° 



for anisotropic modelling IC scattering from PSR B1259- 
63/LS 2883. 

6.2 Modelling the IC gamma-ray spectrum 

The anisotropic IC scattering rate is determined by perform- 
ing the quadruple integration in equation ((6|. This integra- 
tion must be calculated numerically for both the star and 
the circumstellar disc. For the star, the photon distribution 
n p h is assumed to be a blackbody distribution while for the 
circumstellar disc the photon distribution is calculated at 
each point in the orbit using the modified COG method and 
the best fit free-free parameter to the optical and IR data 
(Table Hand Fig.gJ). 

The shape of the electron distribution, n c , depends on 
the location of the electrons. In the pre-shock region it is 
assumed that the electron distribution is mono-energetic 
and only cools through inverse Compton scattering (e.g. 
iKirk et alii 19991 ). The Lorentz factor of the pre-shock pulsar 
wind is typically of the order of j p ~ 10 6 , based on mod- 
elling of the Crab Nebula, whil e the possibility of 7 P ~ 10 4 
electrons has been suggested bv lKhangulvan et all (|2012l ) to 
explain the post-periastron flare observed by Fermi/LAT 
after the December 2010 periastron passage. We have con- 
sidered these two mono-energetic distributions in the AIC 
scattering calculation. 

In the post-shock region the electrons are re-accelerated 
to a power law distribution 

i \ _ \ K c *y~ v 7 m i n < 7 < 7 

max i 

1 elsewhincludedere. 

where K e is a constant and p is the electron index. These 
post-shock electrons can c ool via radiative and adiabati c 
processes. As was shown by I Van Soelen fc Meintiesl (|201ll ). 
if the adiabatical ly cooled electron spectrum considered by 
IKirk et all (|l999l ) is used, the additional target photons as- 
sociated with the IR excess result in an increase in the scat- 
tering rate at GeV energies. This electron distribution has 
an electron index of p — 2.4 and a distribution lying between 
7min = 5.4 x 10 5 and 7 max = 5.4 x 10 7 . We have adopted this 
distribution to model the IC scattering of post-shock elec- 
trons. Additional radiative cooling (e.g. synchrotron emis- 
sion) has not been considered. 

It is not an unreasonable assumption that adiabatic 
cooling dominates near the stand-off shock in PSR B1259- 
63, as the synchrotron radiative cooling time is 
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Table 5. Stellar and disc parameters used to model LS 2883 in 
this study. 




45° disc 



Figure 5. Diagram showing the three disc orientations that we 
considered in this study. The size of the pulsar, Be star/disc and 
orbit are not shown to scale. 



which, if the particles are travelling near the speed of light, 



implies the electrons traverse a distance D a 



367?» be- 



fore cooling. This is comparable to the assumed size of the 
circumstellar disc (i?disc ~ 50R*), therefore near the stand- 
off shock only adiabatic cooling is dominant. 

We considered three orientations of the circumstellar 
disc in this study (Fig. [SJ: the first with the disc plane ori- 
entated at 90° to the orbital plane with the disc normal 
lying in the orbital plane parallel to the semi-major axis; a 
second with the disc plane orientated at 10° to the orbital 
plane, with the sa me rotation; and the last u sing the orien- 
tat ion proposed b y Chernvakova et all l|2006h and modelled 
by lOkazaki et all i^OHl ), with the disc plane orientated at 
45° to the orbital plane and the disc normal rotated 19° to 
the semi-minor axis. Hereafter the orientations will be iden- 
tified as the 90° disc, 10° disc and 45° disc respectively. For 
all cases the circumstellar disc is assumed to have a radius 
R disc = 50i?„. 

For the circumstellar disc, a half-opening angle of 
#disc = 1° was used. For such a small opening angle the 
the angle to the normal, $ in equation can be approxi- 
mated by the polar angle, i.e. $ = 0. The disc temperature 
was taken to be Tdisc = 0.6T* and the disc was considered 
to extend to a distance of i?disc = 50R*. The stellar, orbital 
and disc parameters are summarized in Table [S] 

We considered three electron distributions: a post-shock 
power law distribution with p — 2.4 between 5.4 x 10 5 < 7 < 
5.4 x 10 7 ; a pre-shock mono-energetic distribution centred 
at 7 = 10 4 ; and a pre-shock mono-energetic distribution 
centred at 7 = 10 6 . The IC scattering results from these 
distributions are discussed below. 



Stellar parameters 


M star 


31 M© 


K» 


9.3 R© 


Mp 


1.4 


7* 


33 000 K 


Orbital parameters 


Inclination 


33° 


Eccentricity 


0.87 


Longitude of periastron 


138.65° 


Disc parameters 


T disc (K) 


19 800 


n 


3.055 


Rdisc 


50 -Ratar 


#disc 


1° 


Tilt to orbital plane 90° 


45° 10° 


Tilt to semi-major axis 0° 


19° 0° 



6.2.1 Scattering spectrum - post-shock electrons 

The predicted AIC scattering rate, (dNtatal/dt dei), from the 
post-shocked, adiabatically cooled electron spectrum at pe- 
riastron is shown in Fig. HJ] The solid black line shown in the 
figure represents the scattering rate due to photons origi- 
nating from the (unobscured) star. The dashed lines show 
the contribution due to the circumstellar disc photons, for 
the different orientations considered. Also shown are the 
combined contributions (star+disc) for the different orien- 
tations. The scattering rate has been scaled (normalized) to 
the peak (unobscured) stellar contribution, to illustrate bet- 
ter the ratio between the star and disc contributions. The 
disc contribution is highest at lower energies, since the scat- 
tering occurs in the Thomson limit, while the stellar con- 
tribution peaks at a few x 10 11 eV, due to the scattering 
occurring mainly in the Klein-Nishina limit. 

As is shown in Fig. |6j the 45° and 90° disc contributions 
exceed the stellar contribution below ~ 1 and ~ 50 GeV re- 
spectively. This effectively means that the IC spectrum at 
GeV energies is dominated by scattering from the circum- 
stellar disc for these two orientations. The modelled spec- 
trum shows that the total scattering rate (star+disc contri- 
bution) is significantly greater than the stellar contribution 
alone. 



6.2.2 Light curve - post-shock electrons 

The predicted integrated fluxes around periastron, using 
the parameters above, are shown in Figs. [7] & [5] for the 
Fermi/LAT and H.E.S.S. energy ranges respectively. The 
modelled light curves show the integrated flux around pe- 
riastron, which has been scaled to the maximum (unob- 
scured) stellar contribution. For the Fermi/LAT light curve 
the integration energy range is 0.1 — 100 GeV, while for the 
H.E.S.S. light curves the energy range is 0.380 - 2.321 TeV. 
The H.E.S.S. range was chosen to co rrespond to the results 
presented in lAharonian et all (|2005l ). The lower limit cor- 
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Energy [eV] 

Figure 6. Modelled AIC scattering rate for the post-shock, adia- 
batically cooled electron distribution at periastron. Shown is the 
stellar photon contribution (solid line) and the disc contribution 
(dashed lines) for the different orientations of the circumstellar 
disc. Also shown is the summed contribution from both the star 
and disc (mark lines). The scattering rate has been normalized to 
the maximum stellar scattering rate in the figure. 



responds to their fig. 5, which presents the integrated flux 
above 380 GeV, while the upper limit, 2.321 TeV, is the 
highest reported energy in their fig. 40 

The stellar contribution in Fig. [6] is calculated assum- 
ing the full solid angle of the star is observable by the pul- 
sar/electron cloud (an unobstructed star, applicable to no 
circumstellar disc). However, if the disc is thick enough to 
obscure part of the star, the photon contribution from the 
star will decrease, reducing the stellar contribution to the 
IC scattering. Also shown in Figs. [7] & are the modelled 
gamma-ray light curves if the disc obscures the star, which 
changes the shape and decreases the flux contribution from 
the star around periastron. The gaps in the modified stellar 
light curves correspond to periods when the pulsar is pass- 
ing through the disc. The light curves show that the disc 
contribution is a minimum around the disc crossing epoch 
(r ± 20 days) and a maximum around periastron. 



6.2.3 Scattering spectrum ~ pre-shock electrons 

The modelled normalized AIC scattering spectrum for pre- 
shock electrons with energies 7 = 10 4 and 7 = 10 6 at peri- 
astron are shown in Figs. l9l and 1 101 respectively. Only a 90° 
disc has been considered for the pre-shock electrons. The 
tw o electron energies whic h are considere d, were suggested 
bv lKhangulvan et all (|2012h (7 = 10 4 ) and lKirk et all(ll999f l 
(7 = 10 6 ). The 7 = 10 4 spectrum (Fig. O shows no indica- 
tion of a peak in the scattering rate, usually associated with 
scattering in the Klein-Nishina regime. The highest scatter- 
ing rate occurs at low energies for both the stellar and disc 
contributions, indicating that the scattering occurs in the 
Thomson limit. The predicted scattering spectrum for the 
7 = 10 6 electrons (Fig. 1 10 p shows that the stellar photons 



7 The exact value of the upper limit is given in data released with 
the publication, available from the H.E.S.S. website: www.mpi- 
hd.mpg.de/hfm/HESS/pages/publications/auxiliary/AA4-42-l-htrnl. 




-20 20 
Days from periastron 

Figure 7. Modelled integrated flux for the AIC scattering 
of the post-shock, adiabatically cooled electron distribution at 
Fermi/hAT energies (~ 0.1 — 100 GeV), normalized to the maxi- 
mum stellar contribution. For each disc orientation, i.e. 90° (top), 
45° (middle) and 10° (bottom), the disc contribution is compared 
with the unobscured and the obscured stellar contribution. 




-20 20 

Days from periastron 



Figure 8. Same as Fig.0 but for the H.E.S.S. energy range. 



scatter in the Klein-Nishina regime (implied by the peak in 
the scattering rate at ~ 5 x 10 1 eV), while the disc pho- 
tons scatter in the Thomson regime. The scattering from 
the disc is higher than the stellar contribution for energies 
below ~ 10 10 eV and the inclusion of the disc contribution 
increases the total scattering rate by a factor of ~ 2. 



6.2.4 Lightcurve - pre-shock electrons 

The predicted integrated flux at Fermi/LAT energies around 
periastron are shown for 7 = 10 4 and 7 = 10 6 pre-shock 
electrons in Fig. 111! scaled to the maximum stellar contri- 
bution. Only the 90° disc has been considered, using the 
same orientation as for the post-shock electrons. While the 
7 = 10 4 electrons do not significantly alter the light curve, 
the 7 = 10 6 electrons produce a disc contribution which is 
significant across the period under consideration, resulting 
in a higher flux than the stellar contribution at periastron. 
The disc minimum occurs around the disc crossing epoch 
(~ t ± 20 days) while the maximum occurs around perias- 
tron. 



10 B. van Soelen, P. J. Meintjes, A. Odendaal and L.J. Townsend 





star 




disc 




star+disc — - " 







Energy [eV] 

Figure 9. The normalized scattering spectrum for a mono- 
energetic electron distribution centred at 7 = 10 4 . The figure 
is normalized to the relevant maximum stellar contribution. 




Energy [eV] 

Figure 10. The normalized scattering spectrum for a mono- 
energetic electron distribution centred at 7 = 10 6 . The figure 
is normalized to the relevant maximum stellar contribution. 



7 DISCUSSION 
7.1 Flux increase 

The predicted increase in the flux at periastron due to the 
contribution of photons from the disc is shown in Fie. 1 121 for 
different electron distributions (pre- and post-shock). The 




-40 -20 20 

Days from periastron 

Figure 11. The predicted integrated flux at Fermi/LAT energies 
from the AIC scattering of pre-shocked electrons for 7 = 10 4 (top) 
and 7 = 10 6 (bottom). The figures are normalized to the relevant 
maximum stellar contribution. 
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Figure 12. Ratio of the total flux (star+disc contribution) to the 
flux contribution due to the stellar photons, calculated at peri- 
astron. Indicated on the graph are the adiabatic cooled spectra, 
calculated for the three disc orientations, with and without the 
star being obscured. Also shown are the results for the scattering 
of pre-shock mono-energetic electrons (7 = 10 4 and 7 = 10 6 ). 



predicted increase in th e flux is of the same order a s the 
isotropic approximation (|Van Soelen fc Meintiesl 12011! ) and 
is found to be greater than a factor ~ 2 below GeV energies. 

Since the anisotropic calculation is highly directional, 
photons that arrive over the wide solid angle of the disc do 
not all scatter with the same energy. Unlike the isotropic 
approximation where an average scattering angle was con- 
sidered and the photons were considered to be isotropically 
distributed, the anisotropic model is based on a photon dis- 
tribution which is spread out over the solid angle of the 
disc. This results in a slight decrease in the influence of the 
circumstellar disc in the anisotropic calculation. 

The increase in the AIC spectrum is of the same or- 
der as the isotropic approximation because the IR ex- 
cess predicted in this p aper is higher than considered in 
IVan Soelen fc Meintiesl l|201ll ). This is due to the change 
in the colour excess. A col our excess of E ( B-V) = 1.05 
was previously propo s ed by Ijohnston et alj |l994) . how- 
ever iNegueruela et all (|201ll ) have proposed a newer value 
of E(B-V)= 0.85. Since the reddening effects are lower, less 
correction is applied to the observed flux. Since interstellar 
reddening influences the optical region far more than the IR 
region of the spectrum, this results in a lower optical flux rel- 
ative to the IR flux. This implies that the IR excess is higher 
than previously considered. The disc is therefore brighter be- 
cause of a higher disc temperature and is brighter relative 
to the star. This relative increase in the photon spectrum 
increases the influence of the disc photon contribution. 



7.2 Light curves 

The predicted AIC, gamma-ray light curves (Figs. [7] [8] & 
Hip show that the disc contribution exhibits three proper- 
ties: firstly, the contribution is a maximum around perias- 
tron; secondly, there is a slight increase just before the disc 
crossing epoch and a minimum at the disc crossing; and 
lastly, the rise after the disc crossing is much faster than 
the slow rise just before the disc crossing. This modelled 
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R disc /R * 

Figure 13. Profile of the normalized photon number density, at 
10 12 and 10 14 Hz, from the circumstellar disc at 20 days before 
periastron plotted through the centre of the disc. Note the y-axis 
is scaled logarithmically. 



disc contribution light curve is a combination of the photon 
distribution of the circumstellar disc and the IC scattering 
angle, as discussed below. 

As shown by the profile of the normalized photon num- 
ber density originating from the disc (Fig. I13[) . higher fre- 
quency photons occur mainly near the centre of the circum- 
stellar disc while longer wavelengths remain constant further 
out in the disc. However, the photon number density is much 
higher at optical frequencies, peaking near v ~ 10 14 Hz 
(Fig. [3). These higher frequency photons mainly originate 
near the centre of the disc, close to the surface of the star. 
At the position where the pulsar passes through the disc, 
~ 50i?« from the centre of the disc, the photon density is 
much lower. While the target photons from the disc origi- 
nate from a wide solid angle near the disc crossing epoch, 
the bright central region only extends over a small solid an- 
gle at large polar angles. Near periastron the pulsar is closer 
to the star and the bright region of the disc, providing a 
larger solid angle of brighter photons for IC scattering. This 
results in a much larger disc contribution towards the AIC 
scattering at periastron and not at the disc crossing epoch. 

The scattering from the disc increases rapidly after the 
disc crossing epoch due to more favourable scattering an- 
gles. Due to the system orientation, before the disc crossing, 
all scatterings occur predominantly head-to-tail (0 ac = 0°) 
while after the disc crossing predominantly head-to-head 
(9 SC — 180°) collisions occur. On the other hand, the AIC 
scattering of the stellar photons follows a smooth transition 
which results in a smooth light curve and not the sharp in- 
crease and decrease seen for the disc contribution. 

The combined effects of the photon density and the IC 
scattering angle result in the predicted AIC light curves. 
Additional effects, such as the heating of the circumstellar 
disc due to the pulsar crossing or changes in the disc size 
near periastron, have not been considered. 

It may be possible that shock heating of the circumstel- 
lar disc could locally increase the disc intensity, providing 
additional photons for IC scatteri ng, and significantly in - 
fluence the gamma-ray light curve. iKhangulvan et all (|2012h 
have suggested that the observed Fermi/LAT flare is due 
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Figure 14. Comparison of the predicted model to the 
Fermi/LAT and H.E.S.S. observations. The vertical arrows in- 
dicate upper limits on the Fermi/LAT detections. 

to the shock heating of the disc and the IC scattering of 
pre-shocked pulsar wind with a Lorentz factor of 7 ~ 10 4 . 
However, the model proposes a large opening angle for 
the circumstellar disc and requires the pulsar to leave the 
disc at ~ t + 30 days near the observed Fermi/LAT flare. 
This is in contradiction to the radio observations which de- 
tect the pulsed ra dio emission r + 15 days after periastron 
|Abdo et al.ll201ll ). 

7.3 Consequences of the model 

By considering only the IR excess from the circumstellar 
disc and not considering any of the other consequences of 
changes in the density, the size of the pulsar wind nebula 
with orbital motion or radiative cooling effects, the total 
gamma-ray flux increases at GeV energies by a factor ~ 2. 
While it was thought that the extended size of the disc may 
completely dominate at the disc cross-epoch, it is shown 
by the modified COG method that the IR intensity at the 
distance of the disc crossing is faint and the majority of the 
emission occurs close to the central star. 

Fig. [14] shows the calculated gamma-ray energy spec- 
trum at periastron in comparison to the H.E.S.S. and 
Fermi/LAT observations. The model curve has been been 
linearly scaled to the observations. The predicted gamma- 
ray flux lies orders of magnitudes below the Fermi/LAT 
detection limit and is therefore undetectable with current 
gamma-ray telescopes. At Very High Energies (VHE), while 
the flux is high enough for detection, the predicted variabil- 
ity which would result from a change in the disc IR contribu- 
tion is smaller than the error in the H.E.S.S. measurements 
and therefore would not be detected at a significant level. 
However, Be stars are also known to be variable at optical 
frequencies. If the optical magnitude of the system changes 
significantly, this could result in a higher variability at TeV 
energies. 

7.4 State of the disc post-periastron 

SPH simulations (|Okazaki et all l201ll ; iTakata et ail I2012T I 
of PSR B1259-63/LS 2883 around periastron have suggested 
that the disc could be disrupted by the passage of the pulsar 
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during the periastron passage. As shown by these authors, 
the degree to which the disc is disrupted depends on the 
density of the disc. For a larger density (po ~ 10~ 9 g cm -3 ) 
the disc remains more stable. This could result in changes 
in the observed optical and IR emission from the disc. How- 
ever, observations near periastron on 5 January 2011 with 
VISIR/VLT found a flux that was of the same order as ear- 
lier Band A measurements with the Midcourse Space Exper- 
iment (MSX) (Table |3J. Similarly, spectroscopic observa- 
tions of LS 2883 with the SAAO 1.9m telescope, taken four 
months after periastron show a very prominent Ha emis- 
sion line and no indication of stellar absorption, indicating 
that the circumstellar disc completely dominates the stel- 
lar absorption at the same frequency. These observations, 
of the IR excess and Ha emission line, which took place 
after periastron, show that the Be star is still surrounded 
by a substantial circumstellar disc. The observed equivalent 
width for the Ha line is EW w -50 ± 5 A, very s imilar 
to the measurement reported by Ne gueruela et all (|201 lh 
(EW = —54 ±2 A). This may be an indication that the 
disc remains prominent despite the pulsar passage or, alter- 
natively, may indicate a very rapid recovery. Similarly the 
observations obtained with SALT during May 2012, far away 
from periastron, still show a prominent Ha line. 

7.5 Connection with the Fermi/LAT flare 

The flare detected with Fermi/LAT after the last periastron 
passage was unexpected. We have suggested that the IR ex- 
cess associated with the circumstellar disc can have a large 
influence on the gamma-ray production, increasing the flux 
by a factor ~ 2. However this increase is too low to explain 
the Fermi/LAT flare. As is shown by the COG model, the 
majority of the IR emission is produced close to the stellar 
surface and not near the location of the pulsar-disc crossing. 
This results in the majority of the disc contribution occur- 
ring close to periastron a nd not near the disc cross ing phase. 
It has been suggested bv lKhangulvan et all \20lj ) that the 
pre-shock scattering of 7 = 10 4 electrons could cause the ob- 
served Fermi/LAT flare if the photon density from this disc 
is large enough. These authors have suggested that a lumi- 
nosity of 40 per cent of the stellar luminosity is required, and 
that the additional photons must be originate over a large 
region. However, the observations undertaken during and af- 
ter the previous periastron passage do not seem to suggest 
such a large heating effect. The IR emission from the cir- 
cumstellar disc remains fairly consistent with the 2MASS 
observations, while the Ha observations four months after 
pe riastron is consist e nt wi th earlier observations reported 
by iNegueruela et all ()201 lh . These observations then sug- 
gest that the shock heating which should occur during the 
puls ar-disc crossing will be more localise d aroun d the pulsar. 

IZabalza. Bosch-Ramon Sz Paredesl (|201ll ) discussed 
thermal X-ray emission in gamma-ray binaries as a result 
of the shock heating of the stellar wind. These authors 
estimated the total luminosity of the thermal X-rays as 
some fraction of the total kinetic luminosity of the stellar 
wind. We can make a similar approximation of the total 
amount of energy available for optical/IR emission from 
shock heating by estimating the kinetic luminosity of the 
stellar wind from LS 2883. The total kinetic luminosity can 
be approximated by Lkm ~ Mv 2 , where M is the mass loss 



rate, a nd v is the speed of the wind. Similarly to lKong et al.l 
|201ll ) we can estimate the kinetic luminosity for the polar 
wind to be, 

and for the equatorial wind to be 

rdisc „„ ln3 o ( M \ ( v \ 2 _J 

Lkin ~ 6 - 3xl ° l lO-TM s yr-i j( lO«cms-0 ^ S ' 

The kinetic luminosity of the circumstellar disc can also be 
estimated by considering the disc to be Keplerian, with a 
velocity v = (GAf, /r) 1 ^ 2 . If the density through the cir- 
cumstellar disc is given by equation [T] the mass flow rate 
of the disc through an area AA is given by AM /At ~ pvAA, 
where the density and velocity depend on the distance r 
from the centre of the star. Integrating over the area of 
the disc perpendicular to the direction of the disc motion 
(AA — 2r tan SdiscdA) gives the total kinetic luminosity as 

L kin ~ J P(r)y — - — j 2r tan ^disc dr 

~ 9 x 10 35 erg s _1 . 

Here p has been calculated from equation [1] and the fitted 
COG parameters have been used. Since only a fraction of 
the available kinetic luminosity can be directed into shock 
heating, it seems likely that any increase in the IR/optical 

emission will have a luminosity <C 10 36 erg s" 1 . 

It has also been shown by iKong et alJ l|2012h that 
Doppler boosting of the synchrotron emission could pro- 
duce the required energy level to explain the Fermi/LAT 
flare. The resulting flare occurs because of the preferential 
direction of the shock outflow after the second disc cross- 
ing which results in bulk Doppler boosting of material in 
the cone of the pulsar wind shock front. This suggests that 
the Fermi/LAT flare is the result of the combined effects of 
shock heating and Doppler shiftin g. The SPH simulations 
undertaken bv lTakata et al.l (|2012j 'l show that circumstellar 
disc can be deformed around the pulsar, creating a cavity- 
like structure around the pulsar. Such a shock heated cavity 
may then allow for additional photons, which will allow for 
head-to-head gamma-ray production after the second disc 
crossing. 

Dedicated optical and IR campaigns around future disc 
crossing phases are required to investigate the degree to 
which the circumstellar disc is influenced by the pulsar and 
pulsar wind. 



8 CONCLUSION 

The binary star system PSR B1259-63/LS 2883 has re- 
mained a fascinating source since the detection of the pulsar 
approximately two decades ago. The interaction between the 
Be star and the pulsar produces multi-wavelength emission 
from radio up to TeV gamma-rays, covering 19 decades in 
energy. The optical companion in the system is a Be type 
star, known to possess a circumstellar disc which produces 
an IR flux. This flux is observed as an IR excess which is 
higher than predicted by blackbody and stellar atmosphere 
models. Since the gamma-rays in the system are produced 
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by the up-scattering of photons, the presence of an extended 
and additional photon source could potentially have a sig- 
nificant effect on the gamma-ray production in PSR B1259- 
63/LS 2883. 

We modified the COG method, first proposed by 
IWaters fe Lamersl l| 1984b . which allowed us to predict the 
flux originating from the circumstellar disc at any position 
near the disc. This method has been used to predict the pho- 
ton density as observed from the point-of-view of the pulsar. 
The COG method was constrained using optical and IR data 
obtained from literature and VLT observations. The photon 
density was calculated around periastron and we applied it 
to the problem of AIC scattering. 

The anisotropic scattering considered the change in the 
photon distribution during the orbit, the variation in the 
solid angle subtending the disc, and the direction of scat- 
tering. The variation in the gamma-ray spectrum around 
periastron was calculated for different disc orientations. 

In this study the leptons in the pulsar wind are as- 
sumed to be mono-energetic before the stand-off shock and 
a power-law energy distribution in the post-shock region. 
Both a pre- and post-shock distribution have been consid- 
ered. For PSR B1259-63, it can be shown that adiabatic cool- 
ing dominates near the stand-off shock and, for this reason, 
only adiabatic cooling was considered for the AIC modelling. 

The IC modelling showed that the inclusion of the IR 
excess associated with the circumstellar disc increased the 
GeV gamma-rays by a factor > 2. This relative increase 
in the flux is partly dependent on the choice of the elec- 
tron distribution. For a wide distribution (e.g. 10 4 ^ 7 ^ 
10 7 ), as was shown for the isotropic modelling (fig. 3 in 
IVan Soelen fe Meintiesl l201ll ). or an electron distribution 
centred around 7 ~ 10 4 (Fig. Ill[) . the IC scattering is dom- 
inated by the scattering of stellar photons in the Thomson 
limit. However, this stu dy has shown tha t for the electron 
distribution proposed bv lKirk et al.l l|l999h . which has been 
shown to provide a reasonable fit to the 2004 H.E.S.S. obser- 
vations, the IR excess from the circumstellar disc increases 
the gamma-ray flux by a factor 2 at GeV energies. 

The predicted light curves show (Figs. [7] [8] & QT]) that 
the circumstellar disc is most influential close to periastron 
and, generally, after the disc crossing. This is a result of the 
disc being brightest near the stellar surface and the asso- 
ciated favourable scattering angles after the disc crossing. 
While it was thought the predicted gamma-ray flux may be 
highest near the disc crossing epoch (~ r ± 20 days), the 
low intensity of the disc at ~ 50 stellar radii from the star's 
surface does not provide a sufficient number of photons to 
increase significantly the gamma-ray flux in this region. This 
study has only considered a steady state disc and has not 
considered what effects the pulsar passage may have intro- 
duced on the circumstellar disc. Shock heating of the disc or 
clumping of material may result in regions of higher inten- 
sity which could result in more localised flux increases. 

The additional effects such as changes in the electron 
cooling have not been considered, as these have been ad- 
dressed by other authors. This study has focussed on the 
effect of the IR excess from the circumstellar disc in or- 
der to determine its influence on the production of GeV en- 
ergy gamma-rays. The study shows that the gamma-ray flux 
could be increased by a factor greater than two at GeV ener- 
gies, which is a significant increase. However, the predicted 



gamma-ray flux lies far below that detected by Fermi/LAT 
at the previous periastron passage. A combination of a num- 
ber of effects, including e.g. Doppler boosting, shock heating 
and changes in the shock front are required to explain the 
flare that was observed by Fermi/LAT approximately one 
month after periastron. 
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APPENDIX A: ADDITIONAL CONSTRAINTS 
ON THE COG METHOD 

The specific intensity which orientates along any line-of- 
sight in the disc can be calculated by considering an impact 
vector q relative to the disc. In Fig. [21 the origin of the co- 
ordinate system is at O (the centre of the pulsar), the unit 
vector e* points towards the centre of the disc (and star) and 
ee,d> in the direction of the incoming photon/line-of-sight. 
The value of q is then defined as 

q — fsin£, 

where f is the radial distance in units of R* and £ is defined 
by 

cos£ = e» • ee,^. (Al) 

Any line-of-sight in a direction defined by 9, (j> intercepts 
the disc mid-plane at a position (Figs.[5]&© 

x c i = H tan# cos <j) + L, 
yd — H tan 9 sin (j>, 

where Xd and yd are measured relevant to the disc mid-plane 
such that the origin is in the centre of the star/disc and the 
Xd and yd are measured in the plane of the disc. The height 
of the disc at any point in this coordinate system is 

h d = \Jx 2 +y 2 tan6» disc , 

where 9disc is the opening angle of the disc. For any line-of- 
sight {ee,cj>), we can construct a distance hi which is perpen- 
dicular to the disc mid-plane and intercepts the line-of-sight. 
If hi intercepts the disc mid-plane at any distance Aa; and 
Ay from the point of the disc plane intercept (xd,yd), then 
the length along hi between the disc mid-plane and the line- 
of-sight intercepts is 

hi = ^Ax 2 + Ay 2 tan<9 fc , (A2) 

where 9k = \ — 9. To find the point of the disc intercept Ax 
and Ay can be solved from hi = hd, and therefore 

Ax 2 + Ay 2 tan# fc , = \J ' x 2 + y 2 tan 6> disc 

= \J{x d + Ax) 2 + (y d + Ay) 2 tan di , 

Since in the plane of the disc, Ay — Aa;tan<^> the equation 
above can be solved for Aa; and hence Ay can be calculated. 
The solution gives two values for both Aa; and Ay, which lie 
below and above the disc plane. Special cases occur when 
Aa; or Ay are equal to zero. The correct solution for the 
intercept above and below the plane is determined by the 
direction of the line-of-sight and its orientation to the disc. 
The values of the (xi,yi) and (2:2,2/2) coordinates are then 
given by 

Xl,2 = x d + Axi,2 

yi,2 = j/d + Aj/1,2, 
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and the z coordinates are calculated by 
tan #disc 



z\ = — 



Z2 



■v? 



tan 8a 



Using the solution to the disc intercepts, we then cal- 
culate the angles X\ and X2 in Fig. [2] The two disc inter- 
cepts lie at (x\, j/i, 21) and (3:2,2/2,22) while the centre of 
the electron cloud/pulsar lies at (x — L,y — 0, 2 = //). The 
angles Xi and X2 are the angles between the disc intercepts 
(xi, j/i, 21) & (X2, J/2, 22) and the centre of the electron cloud, 
and are calculated from 

(L,0,H)-(x 1 ,y 1 ,z 1 ) 



cos Xi — 



cos X2 = 



\(L,0,H)\\\\(x 1 ,y 1 ,z 1 )\ 

(L,0,H) ■ (X2,V2,Z2) 



(L, 0,H)\\\\(X2, j/2, 22)|| ■ 
The integration limits on w in equation © are then 



given by 



Wi = 



U12 = 



2 s 



x 2 , 



where £ is defined in equation (IA1[) . 



